Combining dynamical mean-field theory and realistic band structure for V2O3 
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Recent neutron scattering experiments on V2O3 show that the magnetic fluctuations on the metallic 
side of the antiferromagnetic metal-insulator transition are not related to the spin structure of the 
insulator, but rather to the bandstructure-driven spin-density wave phase of the doped system 
V2- y 03. We calculate these magnetic fluctuations starting from a Slater-Koster bandstructure and 
incorporating the correlation effects through the dynamical mean-field theory (DMFT) . Our results 
demonstrate that the magnetic properties of the paramagnetic metallic phase are dominated by the 
Fermi surface topology. On the other hand, the electron-electron interaction drives the paramagnetic 
metal-insulator transition in V2O3. The transition to the antiferromagnetic insulator by virtue of 
orbital ordering is discussed in the framework of the DMFT. 



Vanadium sesquioxide V2O3 is special among the many 
fascinating transition-metal oxides as it served as one of 
the first examples of a material exhibiting an interaction- 
driven Mott-Hubbard metal- insulator transition (MIT). 
The phase diagram of V2O3 (cf. Fig. [I] a) shows four 
different phases. The controlling variable, besides the 
temperature, is the ratio of local Coulomb interaction U 
and bandwidth D as evidenced by the similar response to 
substitutional doping with Cr and Ti or vanadium defi- 
ciency on the one hand and the application of hydrostatic 
pressure on the other hand. At low temperatures, an an- 
tiferromagnetic insulator (AFI) and an antiferromagnetic 
metal (AFM) occur. The spin structures of these mag- 
netically ordered phases differ dramatically [Q . At higher 
temperatures, a first-order paramagnetic metal-insulator 
transition occurs, ending in a second-order critical point. 
The magnetic fluctuations in the paramagnetic phases 
have recently been shown to correspond to the long-range 
order observed in the AFM 0. 

Experimentally, there are indications that the transi- 
tion to the AFM phase is a spin-density wave transi- 
tion, i.e. an electronic instability due to the formation of 
electron-hole pairs because the Fermi surface has nesting 
sheets: Parts of the Fermi-surface are almost parallel, 
providing large phase space of low-energy electron and 
hole excitations at positions k and k + Q in momentum 
space that are connected by the same nesting vector Q. 
We recently demonstrated that this mechanism indeed 
provides an explanation for the properties of the AFM 
phase H). 

The point of departure is the electronic structure of 
V2O3. Early attempts to calculate the bandstructure of 
V2O3 |4j were based on tight-binding calculations in a 
strongly reduced basis set, and are thus at best quali- 
tatively correct. Recently, state-of-the-art density func- 
tional theory (DFT) results were made available ||. To 
incorporate these results we employ the Slater-Koster 
method, which uses the DFT energy eigenvalues at points 
of high symmetry to parameterize tight-binding matrix 
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FIG. 1. Schematic phase diagram of V2O3, a) including 
and b) without the orbital ordering transition to the AFI. 
Cf. text for abbreviations, Ref. || for numerical values. 

elements, which can then be evaluated at any point in 
the Brillouin zone. This reduces the numerical effort to 
obtain reliable bandstructure information. Details of this 
approach are provided in a recent article H. 

Using this bandstructure we are able to calculate the 
properties of the SDW transition || within a weak- 
coupling Hartree-Fock mean-field theory, and find ex- 
cellent agreement with the experimental situation 
The value of the Fermi surface nesting vector responsi- 
ble for the spin-density wave instability is found to be 
Q = 1.69c*, the transition temperature to be 10 K at 
an effective Hubbard interaction U c ~ 0.5 eV, in good 
agreement with optical experiments |t]]. The success of 
the weak-coupling theory shows that the spin-density 
wave is Fermi-surface-driven, and is not a frustration 
phenomenon of the antiferromagnetic insulator as sug- 
gested previously . The occurrence of the paramag- 
netic MIT on the other hand indicates that correlations 
are important in V2O3. 

The dynamical mean-field theory (DMFT) employed 
here is based on the insight that in the limit of high spa- 
tial dimensions or, better, for high lattice co-ordination 
spatial fluctuations become irrelevant, while dynami- 
cal quantum fluctuations remain important as shown 
by Vollhardt and Metzner |iot| . This procedure allows 
to chose one lattice site i = and integrate out all 
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other degrees of freedom, leading to a local effective ac- 
tion S ff[£/oL which is a functional of an effective field 
Qo(t — t'). The effective field describes the interaction of 
the local problem at i = with a bath made up of the 
rest of the lattice. An electron at i — can go into the 
bath at (imaginary) time r and return to the local site 
at t'. Excellent reviews of the DMFT include Ref. (ll). 

The solution of the local problem remains non-trivial. 
Here, we use the so-called modified iterated perturbation 
theory (MIPT), which becomes exact in both the weak 
coupling and atomic limit and at low and high frequencies 
|L2| . The self-energy T,(iuj n ) in MIPT is obtained through 
the Ansatz 
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where T,( 2 '(iu) n ) is the second-order perturbation theory 
result, 
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Details on obtaining a and b to fit various limiting case 
are given in Ref. The lattice enters through a self- 
consistency condition: The effective field of the local 
problem is identified with the local Green's function of 
the lattice problem, and the full lattice Green's function 
G(k, iu> n ) is calculated from the local self-energy. The 
self-consistency equation is 



E(iw n ) = Q Q (iu n ) - G 1 {iuj n ) 



and 



G{iw n ) = ^Gk(iwri) = — 
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The momentum summation brings in the density of states 
of the actual lattice. Before employing the MIPT, we 
point out that Rozenberg et al. |J applied the DMFT 
on a Bethe lattice to get a fairly complete picture of 
the interaction-driven metal-insulator transition, where 
they introduce phenomenological frustration to suppress 
a Slater-type transition. A first central point here is that 
we show that their results regarding the paramagnetic 
MIT essentially carry over to the case of V2O3 including 
the realistic bandstructure as presented in Ref. || . 

In Fig. H we show the development of the spectral 
function A(u>) = — 1/ttG(uj + i0 + ) as the interaction 
U is tuned through the MIT. It is obtained after solv- 
ing the DMFT in the MIPT approximation and Pade- 
transforming the results to the real frequency axis. The 
transfer of spectral weight to higher frequencies and the 
eventual formation of upper and lower Hubbard bands is 
clearly seen. As long as the system stays metallic A(0) 
does not change, as required for a local self-energy. The 
narrowing of the quasi-particle peak at uj = can be 
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FIG. 2. Development of the spectral function as the inter- 
action U is tuned through the MIT. 
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FIG. 3. Divergence of the effective mass m* at the MIT. 
Dotted line is an extrapolation to T — 0. 

described as a divergence of the effective mass m* or 
reduction of the quasiparticle weight Z, given by 

£-3- 1 -5; ,toI! <" + " f >L-- < 5 > 

The result is shown in Fig. |[ The finite temperature 
curves extrapolate to a linear behavior at T = 0. 

Rozenberg et al. found a coexistence regime for metal- 
lic and insulating solutions of the DMFT equations. We 
confirmed this finding for the realistic bandstructure of 
V2O3. Fig. [| shows the stability regimes, the metallic 
solution is stable for U < U C 2 , the insulating solution for 
U > U c \. The physical phase transition occurs at the 
crossing of the free energy surfaces of these curves, in- 
dicated by the dashed line in Fig. |]. The dotted line is 
an extrapolation of U C 2 to T — based on Fig. ||. The 
values of the transition temperatures and the effective 
interaction agree well with experimental values H. 

While the consideration of the realistic bandstructure 
did only provide quantitative corrections of the paramag- 
netic MIT, it is essential for obtaining even a qualitative 
picture of the magnetic properties. Rozenberg et al. ex- 
tended their study to include frustrated magnetism in 
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FIG. 4. Stability regime of metallic and insulating solutions 
to the DMFT equations. Dashed line indicates physical phase 
transition. 

the Bethe lattice, but this approach is not able to explain 
the different spin structures of the AFI and AFM and the 
absence of any AFI-type fluctuation in the PM in close 
proximity to the AFI g]. 

Since the magnetic structure of the AFM phase and the 
corresponding fluctuations in the paramagnetic phases 
are experimentally known to be incommensurate, the 
only way to study them is through the dynamic spin 
susceptibility, as no precursors of a symmetry breaking 
phase can be seen on the one-particle level as d — > oo 
|l3|] . The DMFT is used here as an approximation to 
calculate the momentum dependent susceptibility; it is 
determined by the lattice Green's functions, Eq. (|4|), as 

X°(q,w) =-5^G(k,i/)G(k + q,a; + i/). (6) 

k,i/ 

The full susceptibility x is obtained from \° after the ver- 
tex functions are taken into account |]iT|] , which become 
momentum independent as d — > oo. 

The dynamic susceptibility is connected with the 
neutron scattering intensity through the fluctuation- 
dissipation theorem, which allows the comparison with 
the measurements of Bao et al. [plJTifl . We present three 
characteristic points in the PM phase: A is chosen close 
to the paramagnetic MIT and close to the AFI, B and 
C are at different temperatures above the AFM low- 
temperature phase. In Fig. ^|we show the momentum de- 
pendence of the susceptibility for a given small frequency 
lu. At all temperatures the maximum in the susceptibility 
is at the wave vector Q corresponding to the AFM spin- 
density wave structure. Remarkably, at A no evidence of 
the close AFI phase is seen, in agreement with neutron 
scattering results. The comparison of B and C shows the 
increase of the intensity at Q as the spin-density wave 
transition is approached. 

In Fig. @ the frequency dependence of %(q = Q, to) is 
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FIG. 5. Dynamic susceptibility at different points in the 
phase diagram (cf. text). The wave vector at the maximum 
of x(q, corresponds to the long-range order observed in the 
AFM. 
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FIG. 6. Frequency dependence of the dynamic susceptibil- 
ity at the same parameter values as in Fig. 

shown, scaled approximately on a single curve. This is 
to be compared with the results of Bao et al. |l4| . Bao 
et al. fit their results with the SCR theory for itinerant 
antiferromagnetism Jl5| . Our curve has spectral weight 
transferred to higher frequencies compared to the SCR 
due to the strong correlations, and does improve on the 
fit to the experimental result. 

We demonstrate that the combination of the DMFT 
with the realistic bandstructure provides a natural de- 
scription of the paramagnetic metal-insulator transition, 
the spin-density wave transition, and the magnetic fluc- 
tuations of the paramagnetic phases. This is summarized 
in Fig. |l| b). In particular, this parameter-free calculation 
gives good quantitative agreement with experimental re- 
sults, both for the spin-density wave § and the MIT. The 
absence of any AFI-type fluctuations outside the AFI as 
observed in neutron scattering experiments is due to the 
strongly first order transition to the AFI. 

Let us now turn to the AFI. Castellani et al. sug- 
gested several years ago, that the essential mechanism 
of the AFI is orbital ordering jl7|. In V2O3, the vana- 
dium e g orbitals form a doubly degenerate, quarter-filled 
band. The degeneracy is lifted in second-order perturba- 
tion theory, leading to a favoring of ferromagnetic spin 
and "antiferromagnetic orbital" coupling, i.e. a staggered 
occupation of the two degenerate orbitals. 
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FIG. 7. Proposed cluster of vanadium sites to study orbital 
ordering in the AFI of V2O3. The split circles denote the 
occupation of the two degenerate orbitals at each site | p^ , p7| . 

Castellani et al. were able to apply this idea to V2O3 and 
demonstrate that the observed spin structure of the AFI 
can be explained in this approach. It is thus necessary 
to go beyond the effective one-band model studied so far 

While this complements our picture of V2O3, the ap- 
proach of Castellani has two shortcomings from todays 
point of view: The authors had to rely on the band- 
structure calculations available at the time f| and strong 
correlations have not been taken into account, even in a 
regime where the paramagnetic MIT is known to occur. 
These corrections are not expected to be essential, yet 
for a complete picture of V2O3 it is desirable to include 
them in the calculation. To that end, the DMFT equa- 
tions have to be generalized to a cluster of vanadium 
sites as shown in Fig. 0. This approach is similar to 
introducing two sublattices to study simple AB-type an- 
tiferromagnetism. Here, four neighbor sites of i = have 
to be considered. The Green's functions for these sites 
are related by 



s-iO _ s~iA _ r-iB _ (~<C_ _ (~<D 
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where v is an orbital index and we suppressed the ar- 
gument iuj n . The orbitally dependent DMFT equations 
have been derived by Momoi and Kubo and by Held and 
Vollhardt (l8) for the doubly-degenerate Hubbard model. 
The hopping matrix elements t%j and the eigenvectors 
v4.„ Q (k), where a is a band index, can be obtained from 
the bandstructure calculation. One has 



= E eIk ' (Ri " Rj) E A -( k ) £ k^(k)- 



(8) 



To solve the coupled problem is straightforward in prin- 
ciple but has not yet been performed for V2O3. 

In conclusion, the various phase transitions of V2O3 
are consistently explained once the physics of strong cor- 
relation and the realistic bandstructure are combined. 
Apart from the calculation of orbital ordering within 
DMFT, some other open questions remain. A system- 
atic 1 / d-expansion would complement the assumption of 



high spatial dimensions in the treatment of the strong- 
correlation problem. The importance of disorder, which 
is certainly present at least in the doped V2O3 com- 
pounds, has not been studied. Dobrosavlevic and Kotliar 
]l9[ provided a promising extension of the DMFT to the 
dirty limit. And, finally, the absence of superconductiv- 
ity in V2O3 is still a challenging problem. The solution 
to this last problem promises new insights into the basic 
mechanism of high-temperature superconductivity in the 
cuprates. One possible scenario is that the large local mo- 
ments known to occur in metallic V2- a C>3 Q destroy the 
superconductivity, so that it would be interesting to look 
for superconductivity in completely stoichiometric V2O3, 
where the metallic phase has been stabilized through ap- 
plication of hydrostatic pressure. 
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